Coulomb blockade double-dot Aharonov-Bohm interferometer: giant fluctuations 
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Electron transport through two parallel quantum dots is a kind of solid-state realization of double- 
path interference. We demonstrate that the inter-dot Coulomb correlation and quantum coherence 
would result in strong current fluctuations with a divergent Fano factor at zero frequency. We also 
provide physical interpretation for this surprising result, which displays its generic feature and allows 
us to recover this phenomenon in more complicated systems. 
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Introduction. — As an analogue of Young's double-slit 
interference electron interfering through mesoscopic 
systems, e.g., a ring-like Aharonov-Bohm (AB) interfer- 
ometer with a quantum dot in one of the interfering 
paths, is of interest for many fundamental reasons 
The AB oscillation of conductance has been observed in 
both closed [1] and open geometry , together with ele- 
gant theoretical analysis [5| . Recently, further study was 
carried out for the closed-loop setup, with particular fo- 
cus on the multiple-reflection induced inefficient "which 
path" information by a nearby charge detector Q . 

Going beyond the mere quantum interference, incor- 
poration of Coulomb correlation between the two paths 
should be of great interest. This can be realized by 
transport through parallel double dots (DD) in Coulomb 
blockade regime. For such DD setup, existing studies in- 
clude the cotunneling interference [3, H, [I3| , and two- 
loops (two fluxes) interference with the two dots as an 
artificial molecule [111, [l^]- Remarkably, super-Poisson 
noise and giant Fano Factor were predicted in this sys - 
tem, as generated by the Coulomb correlations [§. [Tsl. Il4| . 

It was very recently found that the Coulomb block- 
ade in parallel dots pierced by magnetic flux $ com- 
pletely blocks the resonant current for any value of $ 
except for integer multiples of the flux quantum $0- It 
was shown there that this effects in a quantum analogue 
of self-trapping phenomenon in non-linear systems. In 
the present paper we concentrate on Coulomb blockade 
effects in parallel dots, where dephasing and lossy chan- 
nels are included. In particular we concentrate on the 
shot-noise spectrum. We demonstrate that in the ab- 
sence of dephasing and lossy channels this quantity di- 
verges at zero frequency. The most important result of 
our analysis is an explanation of this phenomenon us- 
ing symmetry arguments. This explanation displays a 
new way for a simple treatment of complicated Coulomb 
blockade effects in the presence of quantum interference. 

Model. — Consider double dots (DD) connected in par- 
allel to two leads. For simplicity we assume that in each 
of the dots there is only one level, Ei and E2, involved 
in the transport. Also, we neglect the spin degrees of 
freedom. In case of strong Coulomb blockade, the effect 




FIG. 1: (Color online) Schematic setup of a double-dot 
Aharonov-Bohm interferometer. To address dephasing and 
electron-loss effects, a nearby quantum-point-contact (QPC) 
detector and lossy channels (with strength 7i{2)) are intro- 
duced. 

of spin can be easily restored by doubling the tunnel- 
ing rates of each QD with the left lead. The system is 
described by the following Hamiltonian, 

H^Ho + Ht+ + Ud\did\d2 . (1) 

IJ.= 1,2 

Here the first term, Hq = Y.k[^kLa\La.kL + EkRa\^akR], 
describes the leads and Ht describes their coupling to the 
dots, 

Ht = Y^ {t^LdlakL + tf.Ralndf,'^ + H.c. , (2) 

where /i = 1 , 2 and a[.^ and a are the creation opera- 
tors for the electrons in the leads while d\ 2 are the cre- 
ation operators for the DD. The last term in Eq. ^ de- 
scribes the interdot repulsion. We assume that there is no 
tunnel coupling between the dots and that the couplings 
of the dots to the leads, t;jL(i^), are independent of energy. 
In the absence of a magnetic field one can always choose 
the gauge in such a way that all couplings are real. In 
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the presence of a magnetic flux <I>, however, the tunnehng 
amplitudes between the dots and the leads are in general 
complex. We write ipL(fl) = ipL(i?)e**"^<«' , where t^L(R) 
is the coupling without the magnetic field. The phases 
are constrained to satisfy (piL + <Pir — 4>2L — 4'2R. — 4>, 
where </> = 27r$/$o- 

To account for dephasing effect, we introduce a "which 
path" measurement by a nearby point contact (PC) de- 
tector with a model description as in Ref. [l^. To 
make contact with conventional double-slit interferome- 
ter, we also introduce electron lossy channels. Slightly 
differing from Ref. 's', instead of the semi-infinite tight 
binding chain introduced there, we model the lossy chan- 
nels by attaching each dot with an electronic side reser- 
voir, which is particularly suited in the master equation 
approach. The side-reservoir model was originally pro- 
posed by Biittiker in dealing with phase-breaking effect 
|16j . i.e., electron would lose phase information after en- 
tering the reservoir first, then returning back from it. 
But here, we assume that the reservoir's Fermi level is 
much lower than the dot energy. As a result, electron 
only enters the reservoir unidirectionally, never coming 
back. 

Formalism. — The transport properties of the above 
described system, both current and fluctuations, can 
be conveniently studied by the number-resolved mas- 
ter equation [l^, [H, The central quantity of 



this approach is the number-conditioned state, 
of the double dots, where n is the electron number 
passed through the junction between the DD and an 
assigned lead where number counting is done. Very 
usefully, p^^\t) is related to the electron- number dis- 
tribution function, in terms of P{n,t) = Tr[p^")(t)], 
where the trace is over the DD states. From P{n, t) 
the current and its fluctuations can be readily ana- 
lyzed. For current, it simply reads I{t) = ed{n(t)) /dt, 
where {n(t)) = X^n "'^("■i ^'^^ current fluctua- 
tions, we employ the MacDonald's formula, S'(w) — 
2uj dt sin ujt-^[{n'^{t)) — (It)'^], to calculate the noise 
spectrum. Here, {n'^{t)) = Y^^n^Pin^t), and / is the 
stationary current. In practice, instead of directly solv- 
ing P(n, t), the reduced quantity {11? {t)) can be obtained 
more easily by constructing its equation of motion, based 
on the "n" -resolved master equation (Tsl. [Toj. 

Under inter-dot Coulomb blockade, i.e., the DD can 
be simultaneously occupied at most by one electron, the 
Hilbert space of the DD state is reduced to |0) = |00), 
1) = |10), and |2} = |01), where |10) means the upper 
dot occupied and the lower dot unoccupied, and other 
states have similar interpretations. Following Ref. [l9l . 
the "n" -resolved master equation in this basis can be 
straightforwardly carried out as 



(3b) 



(3c) 



Pi2 = e 



^RP22 - 2 (7<i + 27 + 2i A + 2rR) p'{^> (3d) 



P2 



In the above equations, A = i?i — i?2 is the offset of 
the dot levels. Tl(^r) = 27ri:)i(j^)|ti(i^)p, and 71(2) = 
27rZ?i(2)|ti(2)P, are the respective rates for the couplings 
to the left and right leads, as well as to the side reservoirs. 
Dl{r) and Di{2) are the density of states of the leads and 
reservoirs, while tL(i?) and ti(2) are the respective tunnel- 
ing amplitudes. Note that in actual calculation presented 
in this paper we replaced Fl with 2Fl (c.f. [l3,[l^). In 
this work we assume that 71 =72 =7- Finally, 7^ 



characterizes dephasing between the two dots, resulting 
for instance from the "which path" measurement by the 
point contact. 

Note that the master equations (3a)-(3e) include the 
off-diagonal density-matrix elements, so that they ex- 
plicitly display their quantum mechanical nature. How- 
ever these equations can be derived from the many-body 
Scrodinger equation only in the inflnite bias limit in the 
presence of Coulomb blockade [13, [H, [l^ • 
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FIG. 2: (Color online) (A) Current switch and dephasing ef- 
fect under closed geometry (7 — 0) for aligned DD levels 
(A = 0). (B) Phase shift for A 7^ and electron loss ef- 
fect. With increasing the lossy strength 7, the conventional 
double-slit interference pattern is recovered. 



Current. — First, we consider the case without electron 
loss, i.e., 7 = 0. Simple expression for the steady-state 
current is extractable: 



2(7d + 2FR)(I-cos(/)) - 4Asin0 1 
7d(7d + 2FR)+4A2 + To 



(4) 



where Iq = 4rLFR/ (4Fl + Fr), is the current in the ab- 
sence of magnetic flux. However, in the following we will 
use the current of transport through a Coulomb-blockade 
single dot, h = 2FLFR,/(2rL + Fr), to scale the double- 
dot current, in order to highlight the interference fea- 
tures. 

For the limiting case 7d=0, i.e., no dephasing between 
the two dots, from Eq. (g]) we have / = IoA'^/{A'^ + 
Iq [^b. ( 1 — cos 0) — A sin 0] } . Then a novel switching effect 
follows this result: as A ^ 0, / = Iq for (p = 2im, while 
7 = for any deviation of (j) from these values. This re- 
markable behavior can be explained by defining new basis 
states of the DD, (ijj|0) — > rf[,|0}, chosen such that d\\Q) 
is not coupled to the right reservoir, i.e., t2R ^ t2R = 0, 
then the current would flow only through the state d\ 1 0) . 
This can be realized by the unitary transformation [15| 



dj 
d2 



1 

J7 



tlR 



2R 



2R ''IR. 



(5) 



with TV = (^H + ^fi)^^^' which indeed results in t2R. — 0. 
Also, the coupling of rf||0) to the left lead reads 

hd'P) = -e'^^'^-^'^KhLhR e'^ - t2LiiR)/M. (6) 

It follows from this expression that t2L = for ip — 2mr 
provided that tiL/hL = hR/hR, or for (j) = {2n -\- if 



tiL/t2L = —tiR/t2R. Obviously, for noninteracting DD, 
d^|0) has no contribution to current, while d\\0) carries 
a magnetic- flux modulated current. In the case of inter- 
dot Coulomb blockade, however, whether the coupling of 
dj|0) to the left lead is zero becomes of crucial impor- 
tance. If t2L 7^ 0, then the state d||0), carrying the cur- 
rent, will be blocked by the inter-dot Coulomb repulsion. 
As a result, the total current vanishes. However, if the 
state (i||0) is decoupled from both leads, it remains unoc- 
cupied, so that the current can flow through the state 
d||0). As shown above, this takes place precisely for 
tiL/t2L — ±fii?/i2J?- If this condition is not fulfilled, 
the current is always zero, even for (f) — 2Tm. 

As we demonstrated above, the switching effect be- 
comes very transparent in the particular basis of the DD 
states. Still, it is very surprising how such a basis emerges 
dynamically? Indeed, an electron from the left lead can 
enter the DD system in any of SU(2) equivalent superpo- 
sitions of its states. Therefore there exists a probability 
for each electron to enter the DD in the superposition 
that eliminates one of the links with the right lead. When 
it happens, the electron would be trapped in this state. 
Even if the probability of this event for one electron is 
very small, the total number of electrons passing through 
the DD goes to infinity for t ^ 00. Therefore the trap- 
ping event is always realized for large enough time. In 
the presence of Coulomb blockade this would lead to the 
switching effect, as explained above. 

In the presence of dephasing, which is modelled by a 
which-path detection in this work, the switching effect 
discussed above will be smeared out, as shown in Fig. 
2(A). In appearance, the resultant interference pattern 
resembles the usual one of the double-slit interferometer. 
However, both qualitatively and quantitatively, there ex- 
ists remarkable differences, e.g., the unchanged current 
at = 27rn, which is also the fully dephased current. 

As the DD levels deviate from alignment, i.e., A 7^ 0, 
the current switching phenomena will also disappear, as 
shown in Fig. 2(B). Similar to dephasing, from Eq. (4), 
we see that the current at </) = 27rn is unaffected by A, 
too. However, for A 7^ 0, this is not the maximal current. 
Accordingly, a phase shift of the interference pattern is 
implied. In more generic sense, this is nothing but the 
breaking of phase locking (20| . for two-terminal trans- 
port which can appear only under finite bias voltage and 
^pically in the presence of electron-electron interactions 

In Fig. 2(B) we display also the effect of electron loss. 
That is, we introduce lossy channels to make the interfer- 
ometer more and more open. As a result, we see that all 
the above distinguished features disappear and the con- 
ventional double-slit interference pattern is restored by 
increasing the lossy strength 7. The basic reason is that, 
as the dots become increasingly open, the side reservoirs 
would reduce the occupation probability on the dots, thus 
make the Coulomb correlation and back-reflection less 
important. 
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Current Fluctuations. — Current fluctuations are usu- 
ally characterized by the zero- frequency shot noise, which 
can be calculated by the particle-number-resolved mas- 
ter equation approach, using the MacDonald's formula 
as sketched in the formalism. Strikingly, for the present 
Coulomb blockade DD interferometer, we find that the 
zero- frequency shot noise can be highly super-Poissonian, 



and can even become divergent as A ^ 0. In the follow- 
ing we first demonstrate this novel result, then show more 
other features of the noise. 

For coherent DD interferometer, analytical result of 
the frequency-dependent noise can be obtained using the 
MacDonald's formula: 



J 



S{u;) 



[(2rL + r^)A2-(2rL 



I 



A2 



,2\2 



+ 21 . 



(7) 



Here we have assumed </> — 27rn. At zero frequency limit, 
the Fano factor reads 



F 



m 
21 



(4r' 



(2rz 



2^2 



(8) 



Strikingly, as A ^ 0, it becomes divergent! Note that 
this divergence is not caused by the average current /, 
but the zero-frequency noise itself. Very interestingly, 
from Eq. ([7]) , we find that the limiting order of A ^ 
and CO ^ 0, would lead to different results. That is, if we 
first make A ^ 0, then uj ^ 0, the result reads 



F = 



(9) 



which is finite and coincides with the Fano factor of 
single-level transport [l^. The limiting order leading 
to Eq. ([9]) implies that we are considering the noise for 
aligned DD levels. In this constructed above, see 

Eq. (O and Fig. 3(a), the two transformed dot-states are 
decoupled to each other, and one of them also decoupled 
to both leads if </> = 2?™. As a result, equivalently, the 
transport is through a single channel, leading to the Fano 
factor Eq. 

However, for A ^ but ^ 0, the situation is subtly 
different. In this case, the two transformed states are 
weakly coupled, with a strength oc A. Thus, the trans- 
porting electron on state d\\0) can occasionally tunnel to 

(JJ|0), which is disconnected to both leads, and its occu- 
pation will block the current until the electron tunnels 
back to d\\0) and arrives at the right lead. Typically, 
this strong bunching behavior, induced by the interplay 
of Coulomb interaction and quantum interference, is well 
characterized by a profound super-Poissonian statistics. 
In Fig. 3(b), the coarse-grained temporal current with 
a telegraphic noise nature is plotted schematically. We 
see that, as A — > 0, the current switching would become 
extremely slow, leading to very long time (~ 1/A) cor- 
relation between the transport electrons. It is right this 
long-time-scale fluctuation, or equivalently, the low fre- 
quency component filtered out from the current, which 
causes divergence of the shot noise as A ^ 0. This is 



d = iTTn 
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FIG. 3: (Color online) Schematic interpretation for the noise 
divergence. Upper panel: Effective coupling of the DD to the 
leads and between the dots, in the representation of trans- 
formed DD states, i.e., |i) = d\\Q) and |2> = d^|0). Lower 
panel: Coarse-grained temporal current, with a telegraphic 
noise nature which causes divergence of the zero-frequency 
noise when A . 



similar, in certain sense, to the well known 1// noise, 
which goes to divergence as / ^ 0. 

It is quite interesting that a similar divergence of the 
noise-spectrum at zero frequency has been found for 
rather broad conditions (but only for inelastic cotunnel- 
ing regime) in the framework of classical master equa- 
tions. The latter neglects the off-diagonal elements of the 
density matrix and assumes weak enough tunneling 
In contrast, our quantum rate equations approach goes 
beyond these assumptions. On the other hand it shows 
that the switching effect and divergency of the noise- 
spectrum can take place only at the large bias voltage 
[15| . Indeed, by applying the unitary transformation 
one can always decouple one of the states from the right 
reservoir. However, one still needs the total occupation 
of this state sAt^oo. Otherwise the Coulomb blockade 
is not complete and so the switching effect. This condi- 
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FIG. 4: (Color online) Contour plot of the Fano factor versus 
the (scaled) magnetic flux (j> and tunnel-coupling asymmetry 
Tr/Fl, for different DD level detuning: A = F (A), 0.5F (B), 
O.IF (C), and O.OIF (D). 

tion can be realized only in the large bias limit, where the 
energy levels of the dots are far from the corresponding 
Fermi energy. 

Magnetic-Flux Dependence. — The previous study was 
restricted to zero magnetic flux so that only the state 
d\\0) is connected to the leads. Now we proceed to 
nonzero magnetic flux and A ^ 0. In this case, TiL{(j)) 
and T2Li4') &re nonzero in general. By tuning the flux 
from = to TT, the effective coupling T2L is switched 
on, while Til switched off. As a result, the strong current 
fluctuation at = is considerably suppressed owing to 
this transition to transport through d\\0) and S^lO) in 
series. In between = and tt, however, we find a lo- 
cal minimum for the Fano factor, with a common value 
given by Eq. being independent of A. Its location in 
(j), however, depends on A. This is because, for different 
A, one can always find a proper (f>, such that d\\0) cou- 
ples to the left lead, both directly and indirectly through 



^210), with an effective coupling strength T^. Remind 
also that, the coupling of d\\0) to the right lead is ^ Tn. 
Accordingly, the Fano factor of Eq. ^ is reached. 

In Fig. 4, we display the Fano factor versus the (scaled) 
magnetic flux (phase difference 0) and the tunnel- 
couphng asymmetry Tji/Tl. Besides the ^-dependence, 
we see that, with the increase of Tj^/Tl, the Fano fac- 
tor is considerably enhanced and becomes highly super- 
Poissonian. Interpretation for this dependence is referred 
to Ref. ^13!, where the concept of effective fast-to-slow 
channels was proposed. 

Finally, not shown in Fig. 4 includes the effects of de- 
phasing and electron loss. It is clear that, for dephased 
(original) dots, we can no longer construct the super- 
position states d\\0) and (i2|0)- Then, the fast-to-slow 
channel induced bunching behavior is not anticipated, 
and the Fano factor is reduced to the Poissonian value. 
For lossy effect, we conclude that, with increasing the 
lossy strength (7), shorter duration time on dots will 
weaken the role of Coulomb interaction and multiple re- 
flections, making the noise characteristics Poissonian, like 
that from usual random emission. 



Note Added. — After the submission of present work to 
the arXiv:0812.0846-eprint, a very recent paper by Ur- 
ban and Konig was caused into our attention [2l| , where 
the enhancement of shot noise and even divergence were 
found in the absence of inter-dot but in the presence of 
intra-dot Coulomb blockade. In that case, the electron 
spin plays an essential role. In our DD Coulomb block- 
ade regime, however, the electron spin is irrelevant to the 
super-Poisson noise and its divergence. 
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